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Complete breakdown of the Debye model of rotational relaxation near the 
isotropic-nematic phase boundary: Effects of intermolecular correlations in 

orientational dynamics 

Prasanth P. Jose, Dwaipayan Chakrabarti and Biman BagchQ 

Solid state and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, Karnataka, India. 

The Debye-Stokes-Einstein (DSE) model of rotational diffusion predicts that the rotational cor- 
relation times n vary as [1(1 + 1)] _ , where / is the rank of the orientational correlation function 
(given in terms of the Legendre polynomial of rank I). One often finds significant deviation from 
this prediction, in either direction. In supercooled molecular liquids where the ratio n /t*j falls con- 
siderably below three (the Debye limit), one usually invokes a jump diffusion model to explain the 
approach of the ratio ri /n to unity. Here we show in a computer simulation study of a standard 
model system for thermotropic liquid crystals that this ratio becomes much less than unity as the 
isotropic-nematic phase boundary is approached from the isotropic side. Simultaneously, the ratio 
T2/77 (where r\ is the shear viscosity of the liquid) becomes much larger than hydrodynamic value 
near the I-N transition. We have also analyzed the break down of the Debye model of rotational 
diffusion in ratios of higher order rotational correlation times. We show that the break down of 
the DSE model is due to the growth of orientational pair correlation and provide a mode coupling 
theory analysis to explain the results. 



I. INTRODUCTION 

The rotational diffusion model of Debye 0, Q was pro- 
posed originally to explain dielectric relaxation of polar 
molecules and to relate the observed relaxation time to 
viscosity by the use of the Stokes-Einstein relation. The 
Debye-Stokes-Einstein (DSE) model is primarily devised 
to model the Brownian motion in orientational degrees of 
freedom. The DSE model provides the following amaz- 
ingly simple expression for the decay of the l-th rank 
orientational correlation function Cf (t) 



with 



where 



C s l {t)=exp(-t/T l ), 



1(1 + l)k B T/C f 



Cr = 6vV 



(1) 



(2) 



(3) 



In the above expression of rotational friction, V is the 
volume of the molecule in question, rj is the viscosity of 
the liquid, ks is Boltzmann constant and T is the tem- 
perature. The above expressions predict the ratio of the 
first and second-rank rotational correlation times, T1/Y2, 
to be equal to three. The higher rank rotational correla- 
tion times also vary following the [1(1 + dependence. 
In the past, different aspects of DSE relationship have 
been tested namely, (a) the shape dependence |jj 0, , 
(b) the viscosity dependence, and (c) the rank depen- 
dence of the orientational relaxation, in simple as well as 
complex liquids. 
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In many cases, the ratio tx/t% is found to fall be- 
low three. In particular, in supercooled liquids, the 
ratio is known to approach unity at low temperatures 
and related issues (like translation-rotation decoupling) 
that have given rise to considerable amount of discussion 

[ISBIfllillllHEQHEini o ne often 

invokes the breakdown of Debye diffusion model, which 
requires small angle Brownian rotational motion for its 
validity, due to the emergence of large scale hopping in- 
volving large angular jumps. The basic idea is that a 
large jump leads to the decay of C ; s (i) of all rank / at the 
same time, so that all of them have similar correlation 
times corresponding to the average waiting time for this 
large jump to occur. 

The case of r\jr% deserves particular attention be- 
cause it is the oft discussed one. The sensitivity of these 
two correlation times (and hence the respective correla- 
tion functions) to intermolecular interactions is different. 
When the rotating molecule has dipolar interactions with 
the surrounding molecules, one expects the first rank ro- 
tation to be more affected than the second rank one. This 
is a manifestation of intermolecular correlation which is 
reflected in the orientational pair correlation functions. 
Thus, if we consider dipolar spheres (such Stockmayer 
liquid), we should expect the ratio t±/t2 to be larger than 
three. On the other hand, if the intermolecular interac- 
tion has up-down symmetry (as in the case of molecules 
with ellipsoid of revolution), the reverse could be true. 
Thus, one may need to include the role of intermolecular 
correlation to generalize the DSE model appropriately. 

One ideal candidate to test the above argument re- 
garding the role of equilibrium pair correlations is a sys- 
tem of model mesogens undergoing the isotropic-nematic 
(I-N) phase transition. The isotropic phase is both posi- 
tionally and orientationally disordered while the nematic 
phase is still positionally disordered but orientationally 
ordered. Earlier experimental and simulations studies 
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m e m m m m i3 have demonstrated that exis- 
tence of the orientational relaxation of nematogens near 
I-N transition have similarities with that observed in the 
supercooled liquids. The orientational order parameter 
is defined by S = (P2 (cos(9))) , where 9 is the angle of the 
orientation of the molecular axis with the director. In the 
isotropic phase, S = for infinitely large systems while 
it is non-zero (around 0.5) in the nematic phase. The im- 
portant point is that orientational correlation undergoes 
rapid increase as the I-N phase boundary is approached 
from the isotropic side. This large growth in correlation 
would provide a testing ground of the rotational diffusion 
model and the role of intermolecular correlations. 

Different molecular models have been used to test the 
DSE in molecular simulations. In a molecular dynam- 
ics simulation study using the Gay-Berne intermolecular 
potential |25|. de Miguel et al. found the ratio t\/t2 to 
be always less than three in the isotropic phase. They 
found that the departure from the Debye limit was pro- 
nounced when density was reduced or temperature was 
increased; this deviation is a manifestation of the inertial 
decay which gives a ratio [I + l/l] for ti/ti+%. Vasan- 
thi et al. in extensive molecular dynamics simulations 
using the same model have studied the aspect ratio de- 
pendence of the DSE relationship [2^|. Recently Jose 
and Bagchi have studied the breakdown of the DSE re- 
lationship near the I-N phase transition in a system of 
the Gay-Berne ellipsoids of revolution \2$- They have 
shown that the relation between the rotational friction 
and viscosity breaks down near the I-N phase transition. 

The motivation of the present work comes partly from 
recent reports of the observed similarity in the orienta- 
tional relaxation between supercooled liquids and liquid 
crystals. H H HJ |H M 

The deviation of the 
rank dependence from the Debye model (Eq.2) has been 
often discussed in the context of supercooled liquids. As 
discussed earlier, this is attributed to the existence of 
large angular jumps at low temperatures. In this work, 
we investigated the rank dependence of the rotational 
diffusion near the isotropic-nematic phase boundary. We 
find that the ratio T1/T2 becomes much less than unity 



as the isotropic-nematic phase boundary is approached 
from the isotropic side. Simultaneously, the ratio T2/77 
(where r\ is the shear viscosity of the liquid) becomes 
much larger than hydrodynamic value near the I-N tran- 
sition. We have also analyzed the breakdown of the De- 
bye model of rotational diffusion in ratios of higher order 
rotational correlation times. Theoretical analysis shows 
that the breakdown of the DSE model can be attributed 
to the growth of orientational pair correlation. We pro- 
vide a mode coupling theory analysis to explain the re- 
sults. Thus, the present analysis seems to suggest that 
one need not always invoke large scale jump diffusion to 
explain the decrease of the ratio T1/T2 from the Debye 
limit. This view raises some interesting questions which 
we address in the Conclusion. 

In the next section, we describe the system and simu- 
lation details. Results of our molecular dynamics simu- 
lation study of a system of ellipsoids of revolution with 
an aspect ratio equal to three along an isotherm and an 
isochore across the I-N phase transition is presented in 
the section Hill These results show that the ratio T1/T2 
can become much less than unity. This section also in- 
cludes results for higher rank orientational time correla- 
tion functions (OTCF). In section llVl we present a the- 
oretical analysis which can explain theoretical aspects of 
these results. Section Ivl provides a summary of our re- 
sults and concluding remarks. 



II. SYSTEM AND SIMULATION DETAILS 

Here we consider a system of molecules with axial sym- 
metry interacting with the Gay-Berne (GB) pair poten- 
tial that has served as a standard model in the simulation 
studies of thermotropic liquid crystals. In the GB pair 
potential |28l |2jj , each molecule is assumed to be an el- 
lipsoid of revolution having a single-site representation in 
terms of the position of its center of mass and a unit 
vector along its principal axis of symmetry. The GB 
interaction between molecules i and j is given by 



U% B (r t] ,e t , ej ) = Uif^.o,. e,)( Ai / 2 - Pij 6 ) 



(4) 



where 



Pij 



&{rij,ej,ej) + op 

^0 



(5) 



Here <jq defines the cross-sectional diameter, is the distance between the centers of mass of molecules i and j, and 
Yij = Yij/rij is a unit vector along the intermolecular separation vector r^ . The molecular shape parameter a and 
the energy parameter e both depend on the unit vectors and ej as well as on f y as given by the following set of 
equations: 



a(rij,ei,ej) = a Q 



1 - 



X f (e; ■ iij + ej ■ f.y) 2 (ej 



h 3 ? 



1 + x(e, • ej) 



l-x(ei-ej-) 



-1/2 



(G) 
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with x—(k 2 — 1)/ (k 2 + 1) and 



where the exponents [i and v are adjustable, and 



e(f ii ,e i ,e j ) = e [ei(e i; e i )] ly [e 2 (fy, e.,, e^)]^ 



eiCe^ej) - [1 - x (e* ■ erf] 



21-1/2 



and 



e 2 (r i ,,-,e i , &,•) = 1 - — 



l + X'(ei-ej) 1 - x'(ej • e,-) 



(7) 

(8) 

(9) 



with X ' = W 1/M - l)/0s' 1/M + !)■ Here k = CT ee /cr s;5 
is the aspect ratio of the molecule with a ee denoting the 
molecular length along the major axis and a ss = <jq, 
and k! = t ss /e ee , where e ss and e ee are the depth of 
the minima of potential for a pair of molecules aligned 
parallel in a side-by-side configuration and end-to-end 
configuration, respectively. It follows that the GB pair 
potential defines a family of potential models, each mem- 
ber of which is characterized by a set of four parameters 
(k,k',(x, v). In the present work, we employ the Gay- 
Berne pair potential with the original and most studied 
parametrization (3, 5, 2, 1) |3Cj . 

All quantities are given in reduced units defined in 
terms of the Gay-Berne potential parameters eo and <Jo- 
length in units of <To, temperature in units of eo/fcs, and 
time in units of (mer^/eo) 1 / 2 , m being the mass of the 
ellipsoids of revolution. We set the mass as well as the 
moment of inertia of the ellipsoids equal to unity. The 
simulation is run in a microcanonical ensemble with the 
system in a cubic box with periodic boundary conditions. 
Further details of the simulation can be found elsewhere 

M 

The system of Gay-Berne ellipsoids of revolution with 
aspect ratio 3 is studied separately along an isotherm 
with density variation and along an isochore with tem- 
perature variation across the I-N transition. As com- 
pared to the density driven transition, the temperature 
driven I-N transition in the present system is known to 
be rather diffuse [3(J • Next we present the results of our 
study. 



III. RESULTS 

A. Density variation along an isotherm 

In figure Q we show the orientational order parame- 
ter variation with density as a system of 576 Gay-Berne 
ellipsoids of revolution transits across the I-N transition 
along an isotherm at temperature T* = 1. The phase 
transition is found to occur over a range of density be- 
tween 0.305 and 0.315. In figures |2Ja) and[2Ib), we show 
the decay of the single-particle orientational correlation 
function for the first eight ranks in the isotropic phase 
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FIG. 1: The evolution of the orientational order parameter 
with density along the isotherm at temperature T* = 1 



and near the I-N phase boundary, respectively. The Zth 
rank single-particle orientational time correlation func- 
tion (OTCF) is defined as 



cat) 



Ei^te(0).ei(0)) : 



(10) 



where ej(0) is the unit vector along the symmetry axis of 
the zth ellipsoid of revolution. The DSE model is found 
to hold good for all ranks of the single-particle orienta- 
tional time correlation functions shown in Fig. GJa) m the 
isotropic phase. As the rank of the correlation function 
increases, the relaxation time decreases. 

Near the I-N phase boundary, as shown in Fig[3 the 
relaxation of the single-particle OTCFs slows down con- 
siderably for all ranks. However, Cf(t) gets affected dif- 
ferently for different I values. The even and the odd 
l-th. Cf(t) behave differently with the appearance of a 
pronounced plateau in Cf(i) and C|(t). Interestingly a 
similar decay behavior has been observed in supercooled 
liquids HE|. 

Figure 2(b) shows that although the initial decay of 
Cf (£) and C|(i) is faster than that of C{{t) and CJ(i), 




FIG. 2: The time evolution of the single-particle orientational 
time correlation functions, whose rank range from 1 to 8, 
shown in a log-log plot at two densities corresponding to (a) 
p*— 0.285; (b) p* =0.315. The curves are arranged in the 
decreasing order of ranks from the left to the right in each 
plot. 



respectively, the decay becomes slower at longer times. 
The even 1-th correlation functions both develop a rather 
long and distinct plateau. The decay of C|(f) is par- 
ticularly revealing because it shows all the four phases of 
decay - the initial Gaussian, followed by the exponential, 
then the crossover to the plateau and the final exponen- 
tial decay. 

In figure |3 we show the evolution of the ratios t±/t2, 
Ti/ts, and n/r4 as increase in density drives the sys- 
tem across the I-N transition. Note that away from the 
I-N phase boundary in the isotropic phase, the ratios re- 
main close to what are predicted by the Debye rotational 
diffusion model. However, as the I-N phase boundary 
is approached from the isotropic side, it is evident that 
the Debye rotational diffusion model breaks down com- 
pletely. 

In figure we show the ratios between the second 
rank orientational correlation time and the higher rank 
(/ = 3,4) orientational correlation times as a function 
of density across the I-N transition. Note that the I-N 
transition affects the second rank OTCF most due to the 
up-down symmetry of the molecular model studied in the 
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FIG. 3: The density variation of the ratios between the 
first-rank orientational correlation time and the second, third, 
and fourth rank orientational correlation times across the I-N 
transition. The circles represent the data for T\jri, the dia- 
monds for T1/V3 and the triangles for T1/T4,. The ratios are 
scaled by the corresponding Debye predictions. 
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FIG. 4: The ratios between the second-rank orientational 
correlation time and the higher rank orientational correlation 
times across the I-N transition. The triangles represent the 
data for T2/T3 and the circles for T2/T4. The ratios are scaled 
by the corresponding Debye predictions. 

present work. Figure 01 shows a cusp like behavior which 
is well-known in the study of equilibrium critical phenom- 
ena of finite sized systems. Its appearance in dynamics 
suggests the existence of large scale fluctuations in the 
orientational order parameter |3lj | . Since the second rank 
orientational correlation function is associated with the 
optical response of the system which undergoes dramatic 
increase near the I-N transition [3lLl3^| . the most affected 
orientational memory function is the second-order one. 

Microscopically this phenomenon may be understood 
qualitatively in terms of the molecular field theory of 
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FIG. 5: The semi-log plot of the ratio ri/rj versus density 
for different I values. The circles show Til 7] for 1 = 1, the 
diamonds for I = 2, the squares for I = 3, and the triangles 
for I = 4. The filled symbols show the ratios for even values 
of/. 



Maier and Saupe [sil l33l | , where the molecule is con- 
fined in an effective field created by its neighbors. This 
effective potential is given by the expression 



A_l 

V*2 



S(3cos 2 (9 l )-l) 



(11) 



where A is constant independent of the temperature, vol- 
ume, and pressure, V is the molecular volume and 9i is 
the angle between molecular axis with a preferred axis. 
This effective potential grows as the order parameter in- 
creases. Note that a it rotation of the molecular axis 
relaxes Cf(t) but not Cf (t). This effect is manifested in 
the higher order orientational correlation functions also. 
However, a random orientation of the smaller angle less 
than 7r/2 is only required for the relaxation of the Cj s (t) 
with I > 3,4, .. etc. Hence the slow down of relaxation at 
these ranks appears as S becomes significantly large. In 
section IV, we present a quantitative theory to describe 
these effects. 

Another important aspect of the DSE which has been 
a subject of intense study in the literature of supercooled 
liquids is the viscosity (rj) dependence of the time con- 
stant (ti) of orientational correlation function. In figure 
03 we present the semi-log plot of this ratio against den- 
sity. The ratio remains a constant in the isotropic phase 
of the liquid crystals, in agreement with hydrodynamic 
prediction. This is also in accordance with the earlier 
simulation j2]j which shows the ratio between the rota- 
tional friction and the viscosity also shows similar behav- 
ior. 

Note that unlike in supercooled liquids, the growth in 
the viscosity of the nematogens near the I-N transition 
is not rapid [27|. Therefore, the ratio between t\ and r\ 
deviates from the DSE prediction only due to the growth 
of 77. It is also evident that the violation of the DSE is 




FIG. 6: The evolution of the ratio between the first-rank and 
the second-rank orientational correlation times with temper- 
ature across the I-N transition (circles). The inclusion of the 
scaling factor ensures that the ratio is equal to unity in the 
Debye limit as shown by the dot-dashed line. On a different 
scale shown on the right is the orientational order parameter 
variation with temperature (squares). 
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FIG. 7: The temperature dependence of the ratio between 
the first-rank and the Ith rank orientational correlation times 
across the I-N transition for 1 = 3, and 4. The inclusion of 
the scaling factor ensures that the ratio is equal to unity in 
the Debye limit as shown by the dot-dashed line. 



found to be different for the odd and the even values of 
I. In figure [31 the ratio between r; and rj for even I values 
grows more dramatically than that for odd I values. 
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FIG. 8: The temperature dependence of the ratio between the 
second-rank and the lih rank orientational correlation times 
across the I-N transition for I — 3, and 4. The inclusion of 
the scaling factor ensures that the ratio is equal to unity in 
the Debye limit as shown by the dot-dashed line. 



the temperature behavior of the ratios / r/ across the 
I-N transition for I = 3, and 4. While only a small devi- 
ation from the Debye behavior is apparent even at high 
temperatures away from the I-N transition, the onset of 
the growth of the orientational parameter marks a sharp 
increase in these ratios. The results, embodied in Figs. EJ- 
|H1 suggest that orientational correlation, that builds up 
across the I-N transition, plays a key role in deviation 
from the Debye behavior of the orientational correlation 
times. 

The contrast between the study along an isotherm and 
that along an isochore reveals the importance of the role 
played by the intermolecular potential in the breakdown 
of the DSE relation. In the study along the isotherm, 
the free volume that is available for the rotation reduces 
thus leading to the formation of the orientational caging. 
In contrast, when temperature is reduced, the attractive 
part of the inter molecular potential becomes more dom- 
inant and results in the formation of the orientational 
caging that arrests the orientational random walk of the 
molecules. The above difference arises because temper- 
ature variation leads only to small changes in density 
because of the dominance of the repulsive part of the 
potential in determining the liquid structure. 



B. Temperature variation along an isochore 

In this work, we have also studied orientational relax- 
ation in a system of 500 Gay-Berne ellipsoids of revolu- 
tion with the use of the same parameterization but along 
an isochore at the density p = 0.32. The drop in temper- 
ature drives the system from the isotropic to the nematic 
phase with the growth in the orientational order. In par- 
ticular, the single-particle second-rank orientational cor- 
relation function decays with such a long time scale near 
the I-N phase boundary that an integral second-rank ori- 
entational correlation time is beyond the scope of the 
present simulation study. The poor data quality at long 
times with the present run length restricts us from having 
a reasonable fit of the long-time decay. In order to have 
an estimate of the second-rank orientational correlation 
time for the sake of comparison of the rank-dependent 
orientational correlation times, we define a correlation 
time t!(T) that is the time taken for the single-particle 
I— th rank orientational correlation function to decay by 
90% at a temperate T. The prime is used to distinguish it 
from the corresponding integral orientational correlation 
time. 

Fig. shows that the ratio t{/tI follows the Debye 
behavior away from the I-N transition, but the onset of 
the rapid growth of the orientational order parameter 
near the I-N transition induces a marked deviation from 
the Debye limit and the ratio falls rapidly. On the other 
hand, the ratios t\/t\ go through maxima on transit from 
the isotropic phase to the nematic phase for both / = 3 
and I = 4 as shown in Fig. \7\ and the maxima corre- 
spond to the temperature below which the orientational 
order parameter is on the rapid rise. Figure |S] illustrates 



IV. THEORETICAL ANALYSIS OF 
ORIENTATION RELAXATION 

Here we present a mode coupling theory (MCT) anal- 
ysis of the above relaxation behaviour. MCT has a long 
and honorable history in describing dynamics near phase 
transitions |34l l35l 36] . Our starting point of the the- 
oretical analysis is the Zwanzig-Mori continued fraction 
representation of the frequency dependent orientational 
time correlation function, Q(z) [H I^IIM HoT] . 



1 



l(l+l)k B T 



z + 



(12) 



/( z +r,(z)) 



where I is the moment of inertia and Ti (z) is the Laplace 
frequency (z) and rank dependent memory function. The 
latter is determined by the torque-torque correlation 
function. In general, it is very difficult to determine 
this correlation function from first principles, but as a 
first approximation, we would combine input from the 
mode coupling theory with that from the time depen- 
dent density functional theory to obtain an expression for 
the memory function, Ti(z), which can be used to under- 
stand the reasons for the breakdown of the DSE model. 
Near the I-N transition the memory function T(z) can be 
written as a sum of two parts 



^bare 



rf ms (z), 



(13) 



where the bare part of the memory function is assumed 
to be rank and frequency independent. This can be de- 
scribed by the two-body (binary) collision model. Note 
that in conjunction with EqH3 T bare leads to the DSE 
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behaviour. The singular part of the memory function 
contains effects of intcrmolecular correlation and is rank 
dependent. It is given by [iH. H^. Eflj . 

rf na (z) = ^fff jH dt e- 2t J™ dk fc 2 £ 4m(m m (k, 

In the above equation, T/ is the rank and frequency de- 
pendent memory function. This is a function of the 
1,1, m component of the wavevector dependent direct cor- 
relation function Cu m (k) (in the inter molecular frame). 
Fi m (k, t) is the /, m component of the orientation depen- 
dent self- intermediate scattering function. The Fi m (h,i) 
is defined in terms of the spherical harmonics as 

F lm (k,t) = (e k '( r «- r (°))y im (ft,o)y im (^))- (15) 

The single particle position and orientation dependent 
T(r — r', t — t', fi, l~i') memory function is related to the 
torque-torque correlation function through the density 
functional theory by following the fluctuation dissipation 
theorem [37],|4(|- Note that T s l m9 (z) contains the integra- 
tion over the wave vector dependence. The slow down of 
the relaxation of single particle orientational correlation 
function is related to the nature of the component of the 
dynamics structure factor. 

It is important to note that the rotational friction de- 
pends on the rank of the orientational correlation func- 
tion and this friction differs from rank to rank because 
of vastly different wave vector dependence of Fi m (k,t), 
for different I, particularly at low wavenumbers. Near 
the I-N transition, due to existence of large wave length 
fluctuations, the k — > component of the F2o(k,t) un- 
dergoes a very slow decay and this is responsible for the 
slow down of the relaxation of the Cf(i). In this limit, 
the expression for Ff m (k, t) is given by [321 

Fi m (k,t) = Si m (k)e~"£nW t , (16) 

where Si m (k) is the orientation dependent structure fac- 
tor and Dr is the rotational diffusion coefficient. The 
Sim{k) is given by the expression 

Si m (k) = (e k < r -^Y lm (n,0)Y lm (n',t)) (17) 

Near the I-N transition, S2o(k) grows as 1/B 2 k 2 (B = 
£k C 'dfc?^ ) )■ The growth of orientational pair cor- 
relation with the approach of the I-N transition is evident 
in figure|Hl Note that the starting from density p*=0.305, 
the orientational pair correlation function becomes non- 
zero even at large intermolecular separations. This is 
reflected in the rapid growth of S f 2o(fc) as k — > near the 
I-N transition. 

Combination of the above factors provides the follow- 
ing simple expression for the frequency dependence of the 
singular part of the memory kernel 

Tr 9 {z) = A/VS (18) 



1 .2 1 1 1 1 1 1 r 




FIG. 9: The (7220 component of the pair correlation function 
versus the pair separation at several densities across the I-N 
transition. The curves staring from the bottom to the top are 
corresponding to densities between p* = 0.285 and p* = 0.315 
on a grid of Sp* = 0.005. 



where A is a numerical constant Note that such 

a power law dependence is absent from all odd ^-th Ti 
but in principle present in all the even £. However, as I 
increases, the decay becomes increasingly inertial and the 
role of intermolecular correlation becomes weak beyond 
I = 4. But for I — 2 and I = 4, the inverse square root 
dependence of the rotational memory function leads to 
a markedly slower, power law decay, as seen from figure 
2(b). 

Physically, this power law is a manifestation of the 
growth of the pseudo-nematic domains near the I-N phase 
boundary. A particle inside this domain feels a localiz- 
ing potential which makes its rotation difficult. However, 
even within such a domain, a rotation of an individual 
particle by 180° (that is, by ir) is possible because of the 
up-down symmetry of the particle. However, in contra- 
diction to the conceptualization prevalent in supercooled 
liquid literature, here such a half-cycle rotation relaxes 
only the odd rank correlation functions, leaving the even 
ranked ones unchanged. In the MCT description, the 
influence of the localizing potential enters through the 
two particle direct correlation function (cu m (k)) and the 
static structure factor S^ m {k). Just as in the Maier- 
Saupe theory, the present mean-field theory description 
can capture the rank dependence of the effective, localiz- 
ing potential. 

The above mode coupling theory analysis is by no 
means complete, but it provides a semi-quantitative ex- 
planation of the observed rank dependence of the orien- 
tational correlation time near the I-N transition, in terms 
of the rapid growth of equilibrium orientational pair cor- 
relation function. 
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V. CONCLUDING REMARKS 

Let us first summarize the main results presented 
here. We study a system of Gay-Berne model meso- 
gens along an isotherm and an isochore separately across 
its isotropic-nematic phase transition to investigate the 
rank dependent single-particle orientational relaxation 
from the perspective of the Debye behavior. Our results 
demonstrate that orientational correlation that starts 
growing near the I-N transition as it is approached from 
the isotropic side induces a marked deviation from the 
Debye behavior. We present a theoretical analysis of 
our results within the framework of the mode-coupling 
theory. This mode coupling theory analysis provides a 
semi-quantitative explanation of the observed rank de- 
pendence of the orientational correlation time near the 
I-N transition. This explanation does not invoke the ex- 
istence of any large scale angular jump motion which 
randomizes and thereby leads to the decay of correlation 
of all ranks with the same rate. Instead, our analysis 
attributes the non-Debye rank dependence to the rapid 
growth of equilibrium orientational pair correlation func- 



tion. 

As already mentioned, an earlier onset of stretching 
in the even rank correlation functions than in the corre- 
sponding odd rank functions has been observed in super- 
cooled liquids as well |9|, [13| ■ This similarity is indeed 
interesting and deserves further study. 

In view of the results presented here, a comparative 
study between dielectric relaxation (which essentially 
measures the £ — 1 correlation function) and light scat- 
tering or fluorescence depolarization or even NMR (all 
of these measure the I — 2 correlation function) should 
be worthwhile. In fact, a detailed theoretical analysis of 
dielectric relaxation near the I-N phase boundary seems 
to be lacking. 
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